Phase diagram of two-flavor quark matter: Gluonic phase at nonzero temperature 
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The phase structure of neutral two-flavor quark matter at nonzero temperature is studied. Our 
analysis is performed within the framework of a gauged Nambu-Jona-Lasinio model and the mean- 
field approximation. We compute the free energy of the gluonic phase (gluonic cylindrical phase 
II) in a self-consistent manner and investigate the phase transition from the gluonic phase to the 
2SC/g2SC/NQ phases. We briefly consider the phase diagram in the plane of coupling strength 
versus temperature and discuss the mixed phase consisting of the normal quark and 2SC phases. 



PACS numbers: 12.38.-t, ll.30.Qc, 26.60. 



I. INTRODUCTION 

The properties of cold and dense quark matter are of great interest in astrophysics and cosmology. In particular, 
at moderate densities of relevance for the interior of compact stars, quark matter is a color superconductor and has 
a rich phase structure with important implications for compact star physics [l|, 0, H, 0, M, H, S B • 

Bulk matter in the interior of compact stars should be color and electrically neutral and be in /3-equilibrium. In 
the two-flavor case, these conditions separate the Fermi momenta of up and down quarks and, as a consequence, 
the ordinary BCS state (2SC) is not always energetically favored over other unconventional states. The possibilities 
include crystalline color superconductivity and gapless color superconductivity (g2SC) [0, EH G3- However, 
the 2SC/g2SC phases suffer from a chromomagnetic instability, indicated by imaginary Meissner masses of some 
gluons [HI . The instability related to gluons of color 4-7 occurs when the ratio of the gap over the chemical 
potential mismatch, A/S/i, decreases below a value Resolving the chromomagnetic instability and clarifying 
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FIG. 1: The phase diagram of electrically neutral two-flavor quark matter in the plane of Ao and T. At T = 0, the g2SC 
phase exists in the window 92 MeV < Ao < 134 MeV and the 2SC window is given by Ao > 134 MeV. The unstable 
region for gluons 4-7 is depicted by the region enclosed by the thick solid line. The g2SC phase and a part of the 2SC phase 
(92 MeV < Ao < 162 MeV) suffer from the chromomagnetic instability at T = 0. The quark chemical potential is taken to be 
/i = 400 MeV. 
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the nature of true ground state of dense quark matter are central issues in the study of color superconductivity 
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12a l2J, 12a |2fl 123, |2a |2a |3a 131|, 13J, [3J, [3J, |3J, |3J, |33] . (For a three-flavor case, see Refs. 

laEaEIfiaSfil-) 

As an example, in Fig. [TJ we plot the phase diagram of neutral two-flavor quark matter in the plane of the 2SC 
gap at Sfi — (Ao) and temperature (T) [3l|, |32|. (The parameter Ao is essentially the diquark coupling strength.) 
In order to obtain this diagram, we employed a gauged Nambu-Jona-Lasinio (NJL) model, which is the very same 
model that we shall use in this paper. Furthermore, we neglected the color chemical potential. The 2SC/g2SC phases 
and the unpaired normal quark (NQ) phase were included in the analysis. The quark chemical potential was taken 
to be fi — 400 MeV, which is a value typical for the cores of compact stars. The region enclosed by the thick solid 
line is unstable (because gluons of color 4-7 have tachyonic Meissner masses there) and, therefore, should be replaced 
by other chromomagnetically stable phases, for instance, gluonic phases [H, HH- (For more detailed discussions of 
the gluonic phases, see Refs. Ha" 5(|.) Note, however, that we did not consider the global structure of a free energy 
in extracting the unstable region, but only the tendency toward the vector condensation (A 6 ) in the 2SC/g2SC 
phases. A self-consistent analysis of the gluonic phases at T = has recently been done by Hashimoto and Miransky 
[36j and they found that the gluonic phase (strictly speaking, the gluonic cylindrical phase II) exists in the window 
65 MeV < Ao < 160 MeV and is energetically more favored than the 2SC/g2SC/NQ phases in this whole window. 

In this paper, we study the gluonic cylindrical phase II at nonzero temperature and revisit the phase diagram shown 
in Fig. [TJ computing the free energy of the gluonic phase in a self-consistent manner. The result would be useful for 
the phase diagram of QCD, and the compact star phenomenology as well. 

II. MODEL 

In order to study the gluonic phase, we use the gauged NJL model with massless up and down quarks: 

C = ij(ip+ A7> + Gd {$ilsee b Ci?) (i>Ci l5 ee b ^) - \F a ^F a ^ , (1) 

where the quark field tp carries flavor = 1, . . . Nf with Nf = 2) and color (a, (3 = 1, . . . N c with N c = 3) indices, 
C is the charge conjugation matrix; {e) lk — e lk and (e b ) a ^ — e bal3 are the antisymmetric tensors in flavor and color 
spaces, respectively. The diquark coupling strength in the scalar (J p = + ) color-antitriplet channel is denoted by 
Gd- The covariant derivative and the field strength tensor are defined as 

= d^-igA a ^T a , (2a) 
F% = d„Ai - d v Al + g t bc A\A c v . (2b) 

To evaluate loop diagrams we use a three-momentum cutoff A = 653.3 MeV throughout this paper. In NJL-type 
models without dynamic gauge fields, one has to introduce color and electric chemical potentials (fig, and /x e ) by hand 
[47l | to ensure color- and electric-charge neutrality. In /3-equilibrated neutral two-flavor quark matter, the elements of 
the diagonal matrix of quark chemical potentials jl are given by 

f-^ur — f^ug 3/-l 7 
Vdr = fJ-dg = M + Sfl, 
Mub = fi~ Sfl- /i 8 , 

fJ-db = ft + S/J, - /i 8 , (3) 



with 



Sfi u s fi e 

M = M - y + y > Sfi = —. (4) 



s-\ P )=[ {s iv *:_ x ), <.--.) 



In Nambu-Gor'kov space, the inverse full quark propagator S 1 (p) is written as 

* + (So) 
with 

(S+y 1 = 7% + (fl~ <W 3 ) 7 ° + trfAp*, (6a) 
(So)- 1 = 7% ~ (P~ <W)7° - S rKT aT , (6b) 
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and 

<f>- = -iee b ~f 5 A, = -iee b ~f 5 A. (7) 

Here t 3 = diag(l,— 1) is a matrix in flavor space. Following the usual convention, we have chosen the diquark 
condensate to point in the third direction in color space. 

For the gluonic cylindrical phase II, B = {gA e z ) is the most relevant condensate, because the chromomagnetic 
instability related to gluons 4-7 corresponds to the tachyonic mode in the direction of S.JU [H, H(| Besides B, we 
have to introduce a color chemical potential fi 3 — (g-Ag) to ensure color neutrality at B ^ 0. Taking into account 
these condensates, the free energy of the gluonic phase in the one-loop approximation is given by 

Q(A, fj, e , fj, 6 , B, fj, 3 ; n,T) 

(4B 2 fi 3 ^B 2 fx 2 a B 2 



2 5 2 2. 9 2 8 ff 2 



12tt 2 V e e 15 

A 2 1^ / d 3 p 



4G D 2^ J (2tt) 



2Tln(l + e-' 3|et ' 1 ) , (8) 



where (3 = 1/T, the e 's are quasi-quark energies and the sum runs over all particle and anti-particle e a 's. Here, we 
added tree- level contributions from gluons (first line on the r.h.s.), 

fl(tree) = y_ jabc jade A eu ^ ^ 

and electrons (second line on the r.h.s.). Note also that the e a 's depend on the vector condensates through the 
covariant derivatives in the quark propagator l[5]). In what follows, we neglect the color chemical potentials and, 
consequently, the tree-level contributions of gluons. We have carefully checked that their effect on the free energy is 
negligible for realistic values of a s ~ 1 (see also Ref. (36|). 

In this work, in order to remove the ultraviolet divergence in the Meissner screening masses we shall use the following 
subtraction 

n R = Q(A, ^B- M , T) - 0(0, 0, B; 0, 0). (10) 

It is known that this free energy subtraction is not adequate to remove the cutoff dependence of the free energy at 
T > 0. In fact, Eq. (fT0|) leads to positive Meissner screening masses in the normal phase at T > [UH^I. In this 
work we do not go into this problem because this unphysical behavior of the Meissner masses is nothing but a cutoff 
artifact and moreover is negligibly small at n — 400 MeV and at the temperatures of interest (20 MeV at most). 

In order to find the neutral gluonic phase, we first solve a set of coupled equations (the gap equation and the 
electrical charge neutrality condition), 

OA Ofl e 

as a function of B and, then, compute the free energy {]r(B). Finally, the minimum of f2 R (B) determines the neutral 
gluonic phase. (In the following Figs. [3] and [4] we plot the free energy evaluated along the solution of the coupled 
equations (fTTj) .) 



III. NUMERICAL RESULTS 



Figure [2] shows A, 5^ and B in the gluonic phase at T = as a function of An. First, let us note that the results of 
Fig. [5] are in good agreement with those shown in Figs. 1, 2, and 3 of Ref. [3a|, where the color chemical potentials 
M3,8 were treated self-consistently. In Fig. [2] one can see that the gluonic phase exists in the window 

66 MeV < A < 162 MeV. (12) 

The gluonic phase is energetically favored over the 2SC/g2SC/NQ phases in this whole window (see also Fig. 5 of Ref. 
[3^|). One also sees that the phase transition between the gluonic phase and the NQ (2SC) phase at An. = 66 MeV 
(162 MeV) is strongly (weakly) of first order. 
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FIG. 2: The gap parameter A (solid line), the chemical potential mismatch Sfi (dotted line) and the gluonic vector condensate 
B (dashed line) versus Ao in the gluonic phase at T = 0. The quark chemical potential is taken to be \x = 400 MeV. 
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FIG. 3: The free energy fl R (B) as a function of B at T = for A = 75 MeV (solid line), A = 85 MeV (dotted line), 
Ao = 100 MeV (dot-dashed line), and Ao = 140 MeV (dashed line). Note that the free energy is measured with respect to the 
2SC/g2SC/NQ phases at B = 0. The results are plotted for fx = 400 MeV. 



Now let us take a closer look at the free energy at T = 0. Figure [3] shows the behavior of Qn(B) measured with 
respect to the 2SC/g2SC/NQ phases at B = 0. The results are plotted for /i = 400 MeV at several values of Ao. 

In the weak coupling regime, 66 MeV < Ao < 92 MeV, the chromomagnetic instability does not exist in the NQ 
phase (see Fig. [1]). We note that the curvature of Qr{B) at B = 0, 



d 2 n R {B) 



''M 



dB 2 



(13) 



5=0 



can be regarded as the Meissner mass squared d 2 Qr/ 8B 2 \b=o in the 2SC/g2SC/NQ phases, since the solutions of 



Eq. CP satisfy A = A + 0(B 2 ) and /i e 



0(B ) for small values of B, where A and ji e denote their values 



We found that m 2 M is indeed zero in the weak coupling regime. In addition, we observed that, 



at B = [24J, |2 

for small values of B, the system is in the ungapped (A = 0) phase and the free energy behaves like Qr ~ 0(B A ) 
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FIG. 4: The temperature dependence of the free energy (measured with respect to the 2SC/g2SC/NQ phases at B — 0) as a 
function of B for A = 75 MeV (a), for A = 85 MeV (b), for A = 100 MeV (c), and for A = 140 MeV (d). The results are 
plotted for /i = 400 MeV and the values of T are given in MeV. 

However, contrary to the result of Fig. [TJ the free energy has a global minimum at B and the gluonic phase 
is energetically favored over the NQ phase. For 92 MeV < Ao < 162 MeV, one finds tachyonic modes at B = 
because the g2SC phase and a part of the 2SC phase suffer from the chromomagnetic instability and, therefore, are 
unstable against the formation of B. Consequently, the gluonic phase is realized in this region, as expected. For 
strong coupling, Ao > 162 MeV, the 2SC phase is chromomagnetically stable in this regime and the free energy has 
a global minimum at B = 0, though it is not plotted in Fig. [3] 

We now turn to the free energy of the gluonic phase at T > 0. Figure 0Ja) display the temperature dependence 
of the free energy for A = 75 MeV. As T grows, the free-energy gain gets reduced, but the change of the vacuum 
expectation value of B is rather small. As a result, we observe a strong first-order transition from the gluonic phase 
to the NQ phase at T ~ 14 MeV. Note that m 2 M remains positive at any value of T, which is consistent with the 
result shown in Fig. [TJ 

In Fig. HJb), the same plot is displayed for Ao = 85 MeV. At low temperature, like in the case of Ao = 75 MeV, 
the gluonic phase is more favored than the chromomagnetically stable NQ phase. At T ~ 9 MeV, m? M turns negative, 
meaning that the stable NQ phase undergoes a phase transition into the unstable g2SC phase (see Fig. [TJ). The 
gluonic phase is energetically favored until the temperature reaches T ~ 20 MeV. Above this temperature, the g2SC 
phase becomes stable and therefore is favored. 

In Figs. IDJc) andUd), we plot the free energy for the cases of Ao = 100 MeV and Ao = 140 MeV, respectively. 
In both cases, the gluonic phase is energetically favored at low temperature. In contrast, at high temperature, the 
global minima of the free energy are realized at B = 0, which means that, as expected from the result of Fig. [JJ the 
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FIG. 5: Schematic phase diagram of neutral two-flavor quark matter at moderate density in Ao-T plane. The thick solid line 
denotes the line of second-order or weakly first-order transitions and strong first-order transitions are indicated by a thick 
dashed line. In the region enclosed by the thick solid and dashed lines, the gluonic phase is energetically more favored than 
the 2SC/g2SC/NQ phases. 

chromomagnetically stable 2SC/g2SC phases are favored. For Ao = 100 MeV, the phase transition from the gluonic 
phase to the g2SC phase takes place at T ~ 21 MeV. In the case of Ao = 140 MeV, the phase transition takes place 
at T ~ 18 MeV. 

Here, we would like to make a comment regarding the order of the phase transitions. As mentioned above, we 
observed the strong first-order transition (gluonic phase <-> NQ phase) at Ao = 75 MeV. On the other hand, for the 
cases of Ao = 85, 100, 140 MeV, the phase transition (gluonic phase <-> 2SC/g2SC phases) is likely to be of second 
order. However, evaluating the free energy self-consistently near the critical temperatures is not easy and hence we do 
not exclude the possibility of weak first-order transitions. Furthermore, it should be also mentioned that, because of 
the cutoff artifact in Eq. (|10p , it might be impossible to distinguish a weak first-order transition from a second-order 
one. 



We studied the gluonic cylindrical phase II at nonzero temperature. Using the gauged NJL model and the one-loop 
approximation, we computed the free energy of the gluonic phase self-consistcntly and investigated the phase structure 
of the gluonic phase. Although we neglected the color chemical potentials, we have checked that, for a s ~ 1, their 
effect on the free energy is negligible. 

In the weak coupling regime, we found that the gluonic phase undergoes a strong first-order transition into the NQ 
phase as it is heated. This is a new aspect of the gluonic phase at T > 0, which is not shown in Fig. [1] On the 
other hand, since the phase transitions from the gluonic phase to the chromomagnetically stable 2SC/g2SC phases 
are of second-order or weakly first-order, we expect that the corresponding critical line shown in Fig. [T](i.e., the right 
branch of the thick solid line) is not drastically altered by the self-consistent analysis. (In other words, the Meissner 
masses squared can be a rough criterion for choosing the energetically favored phase in this regime.) We thus are able 
to make a sketch of a schematic phase diagram of two-flavor quark matter, which is free from the chromomagnetic 
instability related to gluons 4-7 (see Fig. 1 The low-temperature region of the g2SC phase and a part of the 
2SC phase is replaced by the gluonic phase. Furthermore, the gluonic phase wins against a part of the NQ phase and 



It is interesting to note that the neutral single plane-wave Larkin-Ovchinnikov-Fulde-Ferrell state [48l , |49| has a similar phase structure 
as the gluonic phase jS E l3ll . [3^ . [3^1 . 



IV. SUMMARY, CONCLUSIONS, AND OUTLOOK 



A. Phase diagram 




100 150 

Aq [MeV] 



100 



A, [MeV] 



150 



FIG. 6: The free energy of the neutral 2SC/g2SC phase (solid line), the gluonic phase (dotted line), and the mixed phase 
(dashed line) measured with respect to the NQ phase as a function of A for a = 3 MeV/mi (left) and a = 10 MeV/mi 
(right). The three dots on the solid line (Ao = 92, 134, 162 MeV from left to right) denote the edge of the g2SC window with the 
normal phase, the phase transition point between the 2SC and the g2SC phases, and the critical point of the chromomagnetic 
instability. The quark chemical potential is taken to be /j, — 400 MeV. 



enlarges its region. We argue therefore that the gluonic phase which could resolve the chromomagnetic instability 
related to gluons 4-7 is a strong candidate for the ground state of a neutral two-flavor color superconductor in the 
intermediate coupling regime. Alternatives include other types of the gluonic phases [2J, [2f| [3y, [37| , the crystalline 
phases [HI ELS] and the mixed phase It should be mentioned that, at T — 0, the gluonic color-spin locked phase 
is more stable than the gluonic cylindrical phase II in some region of Ao and moreover is free from the chromomagnetic 
instability at moderate densities [H], HH • 

Although we concentrated on the phase diagram in T-Aq plane in this work, it is obviously worthwhile to revisit 
the phase diagram in T-/J, plane. A preliminary study (5(| indicates that currently known phase diagrams 
must be significantly altered. In addition, the critical temperature for the gluonic phase could reach a few tens of 
MeV and, therefore, it is interesting to study astrophysical implications of the gluonic phase, e.g., the quark matter 
equation of state, neutrino emission from compact star cores, and so on. 



B. Gluonic phase versus mixed phase 

Finally we briefly look at a mixed phase consisting of the NQ and the 2SC phases [l|| [2(J HH • For the mixed phase 
to exist, it must satisfy the Gibbs conditions, which are equivalent to chemical and mechanical equilibrium conditions 
between the NQ and the 2SC phases. These conditions end up as follows 

p( N QW e )=p( 2SC W). (14) 

Beside Eq. (|14p two components must have opposite electrical charge densities. Otherwise a globally neutral mixed 
phase could not exist. We solved Eq. (fT4"]) and found that the globally neutral mixed phase exists in the window 

67 MeV < A < 201 MeV, (15) 

where the quark chemical potential was taken to be /i = 400 MeV. 

In order to calculate the free energy of the mixed phase we take account of finite-size effects, i.e., the surface and 
Coulomb energies associated with phase separation. The surface and Coulomb energy densities are given by 

fs = ec _ 2iTa c - m fd{x)x(An e ) 2 rl, (16) 

where a is the surface tension, x is the volume fraction of the rarer phase, An e is the difference of the electric charge 
density between NQ and 2SC phases, and a C m = 1/137. These energy densities also depend on the dimension d 
(d = 1, 2, and 3 correspond to slabs, rods, and droplets configurations, respectively) and r , which denotes the radius 
of the rarer phase. The geometrical factor fd(x) is given by 
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Minimizing the sum of es and ec with respect to tq, we obtain 

es + e c = l (^a cnl d 2 f d (x)x 2 ) 1/3 (Antf'W*. (18) 

The actual value of the surface tension in quark matter is poorly known, in this work we assume d = 3 (droplets 
configuration) and try relatively small surface tension. 

Figure [6] displays the free energy of the 2SC/g2SC phase, the gluonic phase, and the mixed phase. For a very small 
surface tension a = 3 MeV/fm , the mixed phase is the most favored in a wide range of Ao, 103 MeV < Ao < 166 MeV. 
The gluonic phase is energetically more favored than the mixed phase only in the weak coupling regime. Note that 
the value of the surface tension, a = 3 MeV/fm at fi = 400 MeV, is close to that calculated by Reddy and Rupak 
(l9| . For a surface tension a = 10 MeV/fm 2 , there still is a wide window where the mixed phase is more stable than 
the g2SC phase, but the mixed phase is less favored than the gluonic phase. 

It should be mentioned here that, however, we did not take into account the thickness of the boundary layer, which 
has been estimated to be comparable to the value of the Debye screening length in each of the two phases, and 
therefore the results shown in Fig. J6]is not a final conclusion (HH, |SZ| • The effect of charge screening would increase 



the surface energy substantially 58, 5 
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